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We report on infrared, Raman, magnetic susceptibility, and specific heat measurements on 
CdCr2 04 and ZnCr2 04 single crystals. We estimate the nearest-neighbor and next-nearest neighbor 
exchange constants from the magnetic susceptibility and extract the spin-spin correlation functions 
obtained from the magnetic susceptibility and the magnetic contribution to the specific heat. By 
comparing with the frequency shift of the infrared optical phonons above Tn, we derive estimates for 
the spin-phonon coupling constants in these systems. The observation of phonon modes which are 
both Raman and infrared active suggest the loss of inversion symmetry below the Neel temperature 
in CdCr204 in agreement with theoretical predictions by Chern and coworkers [Phys. Rev. B 74, 
060405 (2006)]. In ZnCr204 several new modes appear below Tn, but no phonon modes could be de- 
tected which are both Raman and infrared active indicating the conservation of inversion symmetry 
in the low-temperature phase. 

PACS numbers: 75.40.-s, 75.50.Ee, 78..30.-j 



I. INTRODUCTION 

The antiferromagnetic oxide-spinel systems ^Cr204 
{A = Cd, Mg, Zn) are prototypical examples for highly 
frustrated magnetSfi where the magnetic Cr'^+ ions with 
spin 3 = 3/2 reside on the vertices of corner-sharing tetra- 
hedra forming a pyrochlore lattice as shown in Fig. [TJ 
An antiferromagnetic nearest-neighbor (nn) Heisenberg 
exchange on the pyrochlore lattice leads to inherent frus- 
tration and a multiply degenerate magnetic ground state. 
In real systems this degeneracy will be released at fi- 
nite temperatures by coupling to other degrees of free- 
dom and, indeed, the ^Cr204 systems {A = Cd, Mg, 
Zn) undergo a magnetostructural transition with antifer- 
romagnetic order at 7.8, 12.7, and 12.5 K^i^ although 
their Curie- Weiss temperatures (Qcw) are —71, —346, 
and —390 K, respectively.^^ The theory of the magne- 
toelastic coupling and the magnetostructural transition 
for these systems was worked out by Tchernyshyov and 
coworkers in terms of a spin-driven Jahn- Teller effect)^ 

The structural and magnetic ground state properties 
in the Cr spinels are strongly intertwined and, hence, the 
assignment of the correct symmetry groups is a complex 
task which has not been accomplished yet: The high- 
temperature cubic structure of spinels has space group 
Fd3m, which is often reduced to tetragonal with space 



group /4i /am(i^>^>^>i^ when undergoing a structural or 
magnetostructural transition. Such a uniform tetrago- 
nal distortion was reported for CdCr204 J ^" "' "' However, 
it was suggested from a theoretical analysis of the ob- 
served spiral spin structure that the lattice distortion 
should be chiral with space group /4i22, leading to the 
loss of inversion symmetry (IS)j^ The corresponding dis- 
tortion in ZnCr2 04 has been mostly described in terms 
of the tetragonal /4m2 space group also lacking inver- 
sion symmetry,— d^d^ but orthorhombic distortions have 
also been evoked leading to F222 symmetryi^ (no IS) or 
Fddd symmetry (conserving IS).— 




FIG. 1: (Color online) The pyrochlore structure formed by 
the magnetic Cr ions on the -B-site in AGvzXi (X = S, O) 
spinels. 
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The number of infrared (IR) and Raman active optical 
phonons strongly depends on the exact symmetry and 
occupied Wyckoff positions of the lattice. If inversion 
symmetry is present IR and Raman modes are mutually 
exclusive. Hence, the appearance of resonances which 
are IR and Raman active is a secure indication for the 
loss of inversion symmetry.— While we are not aware 
of any low-temperature Raman studies of these systems, 
the magnetoelastic coupling and the phase transitions 
clearly show up in the far-infrared optical spectra as shifts 
or splittings of the phonon modes . "^'^'^^'^^ For ZnCr204 
the spin-phonon coupling was estimated via two differ- 
ent approaches: The first one is a direct evaluation of 
the phonon splitting in the magnetically ordered phase 
using the spin-driven Jahn- Teller effect and the second in- 
volves a comparison of the frequency shift of the phonons 
with the spin-spin correlation function derived from the 
magnetic contribution to the specific heat.— The first ap- 
proach also has been applied to CdCr2 04 alreadyA^ 

In this study we investigate the low-temperature Ra- 
man and IR optical response for both systems and an- 
alyze in detail the approach to derive the spin-phonon 
coupling by considering both, the spin-spin correlation 
function derived from the magnetic contribution to the 
specific heat, and the one obtained directly from the 
paramagnetic susceptibility. The outHne of the paper 
is as follows: In section IIII Al we describe the suscepti- 
bilities of CdCr204 and ZnCr204 by a Quantum Tetra- 
hedral Meanfield model^fi and obtain the nn and next- 
nearest neighbor (nnn) exchange couplings in these sys- 
tems. Then we extract the magnetic part of the specific 
heat in Sec. IIIIBI After discussing the phonon excita- 
tions observed by IR and Raman measurements with re- 
gard to the proposed lattice symmetries in section BlI C| 
we estimate in Sec. IIIIDI the spin-phonon coupling in 
both compounds by comparing the frequency shift of the 
IR active phonons to the spin-spin correlations obtained 
from the susceptibility and the specific heat data. 



II. EXPERIMENTAL DETAILS 

High-quality platelike single crystals of CdCr204 up to 
7 mm in diameter were obtained by spontaneous crystal- 
lization from Bi203-V205 fiux. ZnCr204 single crystals 
were grown by a chemical transport technique using CI 
as a transport agent at temperatures between 950 and 
900° C. The length of the octahedra edges extents up to 
6 mm. The optical measurements were performed on as 
grown (111) mirror-like surfaces. X-ray powder diffrac- 
tion patterns of crashed single crystals did not reveal 
any impurity phases and correspond to stoichiometric 
polycrystals. The magnetic susceptibility measurements 
were carried out in a commercial superconducting quan- 
tum interference design magnetometer (Quantum Design 
MPMS-5). The heat capacity was measured in a Quan- 
tum Design Physical Properties Measurement System for 
temperatures from 2 K < T < 300 K. The IR refiec- 
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FIG. 2: (Color online) Temperature dependence of the mag- 
netic susceptibility of CdCr2 04 and ZnCr2 04 in a field of 
fioH = 1 T and ^oH = 0.1 T, respectively. The solid lines 
are fits to the experimental data according to Eq. 



tivity was determined using a Bruker Fourier- Transform 
IR spectrometer IFS 113v/S, which was equipped with a 
He-bath and a He-fiow cryostat (5 K< T < 300 K). Ra- 
man scattering experiments were carried out in a quasi- 
backscattering geometry. An Ar/Kr ion laser was used 
for excitation at 514.5 nm (2.41 eV). The laser output 
power was kept below 3 mW on a focus of approximately 
50 fiui of diameter to protect the sample from possi- 
ble heating effects. The scattered fight was collected 
and dispersed by a triple monochromator DILOR XY 
on a liquid-nitrogen-cooled CCD detector in parallel and 
crossed polarization configurations. Temperature depen- 
dencies were obtained in a variable temperature closed 
cycle cryostat (Oxford/Cryomech Optistat, RT-2.8 K). 

The measured refiectivity spectra are analyzed by 
modeling the dielectric function and fitting it to the ex- 
perimental data using the program ReFFit.— We utilize 
a three-parameter model for each phonon mode which is 
a sum over Lorentzian oscillators: 



(1) 



Here, flj denotes the plasma frequency, woj the eigen- 
frequency, and the damping of mode j, respectively. 
Soo is the dielectric constant beyond the phonons in the 
mid-infrared regime and accounts for high-frequency elec- 
tronic polarizations. It is treated as a free parameter in 
the fitting procedure. The optical phonons in CdCr204 
and ZnCr204 were measured from 5 to 300 K in the far- 
infrared and additional measurements were performed in 
the mid-infrared up to 8000 cm~^ for CdCr204 to de- 
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termine the value of Eoo = 4.6 with high accuracy and 
to enable a Kramers-Kronig transformation. Therefor a 
constant extrapolation of the measured reflectivity was 
assumed at low frequencies and a smooth uj~'^ behavior 
at the upper limit of the recorded spectrum. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Magnetic Susceptibility 

Figured shows the magnetic susceptibility of CdCr204 
and ZnCr204 as a function of temperature. On cooling 
both curves exhibit a broad maximum before magnetic 
order sets in at the indicated Neel temperatures. Note 
that the susceptibility below T/v was shown to depend 
strongly on the applied magnetic fieldi^ Our results are 
in good agreement with previous studies i^i^^i 

Garcia- Adeva and Huber derived a Quantum Tetrahe- 
dral Meanfield (TMF) model^S, to calculate the magnetic 
susceptibility of pyrochlore lattices. We use this model 
in the form 



Xtmf{T) 



a ■ Xtet{T) 



l + (3Ji + 12J2)xtet(T)' 



(2) 



which describes the susceptibility per mole of magnetic 
ions. Here, Na is the Avogadro constant, g = 1.97 the g- 
factor,— the Bohr magneton, and Ji and J2 are the nn 
and nnn exchange constants (in units oiks)-, respectively. 
Xtet can be calculated through 



E g{S)S{S + 1){2S + l)e "'^2^"^" 



XtetiT) = 



12T 



s 



(3) 



The sum runs over the total spin values S = 
(0,1,2,3,4,5,6) of the Cr-tetrahedron and g{S) — 
(4, 9, 11, 10, 6, 3, 1) are the corresponding degeneracies.— 
The scaling factor a is related to the effective number of 
Bohr magnetons by a = p1ff/[g'^s{s + 1)]. 

The solid lines in Fig. [2] were obtained by fitting 
Eq. ^ to the experimental data. This procedure yielded 
Ji = 14.7 K, J2 = -4.0 K, and Peff = 4.04 for CdCr204 
and Ji = 33.4 K, J2 = 4.4 K, and Pe// = 3.69 for 
ZnCr204. In literature values for Pe// range from 3.98 
to 4.02 for the Cd spinel^i"^ and between 3.85 and 
3.94 ^B for the Zn compound.— »24 The spin only value 
for a Cr"^+ ion of 3.87/iB corresponds well to the reported 
data and to our findings. The exchange coupling con- 
stants were estimated from the Curie- Weiss temperatures 
to be Ji = 12 K (Ref. ^ for the Cd and Ji = 35 - 45 K 
(Ref. [23} ) for the Zn compound which is also in good 
agreement to our results. In Ref. I25I the TMF model was 
already applied to ZnCr204 and yielded Ji — 39.4 K and 
J2 = 1.76 K. Evidently, Ji is about one order of magni- 
tude stronger than J2 in ZnCr2 0, while in CdCr204 the 
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FIG. 3: (Color online) Upper frame: Temperature de- 
pendence of Cp/T of CdCr2 04 (open circles), CdGa2 04 
(open squares), and the reference compound ZnAl2 04 (open 
rhombi) together with the fitted lattice contributions C(^l /T 
of the specific heat (lines) as described in the text. Upper 
insets: Cp/T as a function of temperature around the AFM 
transition. Lower frame: Cp/T of ZnCr2 04 vs. temperature 
(open triangles) together with a fit of the phonon contribution 
(dashed line). Also shown are ZnGa2 04 (open hexagons) and 
the reference compound. Lower inset: Entropy of the mag- 
netic contribution Cp — C{^^ of the specific heat for both Cr 
spinels. 



antiferromagnetic Ji is already weakened and the nnn ex- 
change J2 is ferromagnetic. This is in agreement with the 
reduction of the direct exchange between nn Cr ions with 
an increase of the lattice parameter when substituting Zn 
by Cd.-i 

We also attempted to analyze our susceptibility data 
with an alternative model proposed by Garcia- Adeva and 
Huber,— >2^>21i^ which, however, did not allow to fit the 
data satisfactorily.— 



B. Specific Heat 

In Fig. [3] we show the temperature dependence of spe- 
cific heat divided by temperature Cp/T for CdCr204, 
ZnCr204 and the non-magnetic compounds CdGa204, 
ZnGa204, and ZnAl2 04. The magneto-structural tran- 
sitions of CdCr204 at Tjv = 7.8 K and ZnCr204 at 
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Tjv = 12.5 K are clearly visible as very sharp peaks, 
indicating first-order phase transformations. 

In order to extract the magnetic contribution to the 
specific heat we followed the procedure used in Ref . [H for 
ZnCr204: We prepared the non-magnetic reference com- 
pounds CdGa204 and ZnGa204 to subtract the lattice 
contribution to the specific heat directly. The obtained 
specific heat data are also displayed in Fig. [3l Obvi- 
ously, the data for the reference compounds do not allow 
for a direct subtraction, and it is not easy to explain 
the discrepancies for ZnGa204 in comparison to Ref. HI. 
Therefore, we chose an alternative approach to extract 
the magnetic contribution. As a non-magnetic reference 
we prepared ZnAl2 04 and fitted the specific heat data 
by modeling the phonon density of state (PDOS). We 
were able to describe the experimentally obtained val- 
ues with one isotropic Debye- (D) and three isotropic 
Einstein-type (i?i,2,3) contributions satisfactorily (com- 
pare rhombi and soHd line in Fig. [3]) with the correspond- 
ing Debye and the Einstein temperatures 9d and 9ei^2 3) 
respectively. To comply with the 21 degrees of freedoms 
per formula unit, we fixed the ratio between the Debye- 
and Einstein-terms to D : Ei : E2 : ^ 1 : 1 : 3 : 2. 
The parameters obtained by a best fit are 9d ^ 330 K, 
0E, = 350 K, Oe^ = 640 K, and 0e, = 1080 K for 
ZnCr204. The derived Einstein temperatures are in good 
agreement with maxima in the PDOS reported by neu- 
tron diffraction studiesi^ 

Turning now to the frustrated systems ZnCr204 
and CdCr204 we modeled the phonon contribution by 
sc aling the Deb ye temperature of our reference with 
Y^tz^aW^^acc) as it is expected for the phonon disper- 
sion of acoustic modes close to the center of the first 
Brillouin zone. mzAO denotes the molar mass of the ref- 
erence and mAco of the Cr spinel, respectively. The 
Einstein terms were adapted in such a way, that the ob- 
tained magnetic contribution Cm — Cp — Cf^* results in 
a magnetic entropy (see lower inset of Fig. [3]) 



Sm{T) = f^d^, 



(4) 



which increases continuously with temperature and ap- 
proaches the expected high-temperature limit value of 
2R ■ ln(2s -1-1) for a spin system with s — 3/2. For 
CdCr204 this procedure yielded 6d — 265 K, 6ei — 
300 K, 0E2 = 570 K, and 0Es = 820 K. In ZnCr204 
we obtained 0d = 290 K, 0e^ = 360 K, 0e, = 600 K, 
and 0E3 = 860 K. 



C. Infrared and Raman modes 

1. Infrared spectroscopy 

In Fig. mja) we show refiectivity spectra for CdCr204 
at 5 K and 15 K, below and above the magnetostructural 
transition at T/v = 7.8 K. The fits using Eq. ^ are 
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FIG. 4: (Color online) (a) Reflectivity spectra of CdCr204 
in the paramagnetic state at 15 K (shifted for clarity) and 
at 5 K, below Tn. The Tiu{2) mode is clearly split at 5 K.— 
(b) Splitting of the Tiu{2) phonon in ZnCr204. (c), (d) Tem- 
perature dependence of the phonon eigenfrequencies of the 
Tiu{2) modes revealing the splitting at Tm in CdCr2 04 and 
ZnCr2 04, respectively. 



shown as soHd lines in the same plot. The four IR active 
phonon modes Ti„(j) (with j = 1, 2, 3, 4) observed above 
T/v (150, 365, 465, and 605 cm~^ at room temperature), 
were labeled according to the triply degenerate modes 
expected from the irreducible representations in the cubic 
Fd3m (#227) symmetry;^ 
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For ZnCr2 04 the refiectivity spectra exhibit a qualita- 
tively similar shape and have been pubHshed earlier.— The 
four reststrahlen bands are located at room temperature 
at 186, 373, 509, and 620 cm~^, respectively. Therefore, 
we focus on the mode Ti„(2) at 8 K and 20 K, below and 
above the Neel temperature Tn = 12.5 K in Fig.IDJb). As 
previously reported^ii^ this mode softens upon cooling 
and splits for both compounds when entering the mag- 
netically ordered phase [see Figs.lDJc) and (d)]. The soHd 
line in Fig. [4][c) is a fit of the increase of the phonon fre- 
quency with decreasing temperature using 



ujj{T) = woj 



1 - 



exp(ft,c LUav/ksT) ~ 1 



(5) 



to include anharmonic effects.— Here, Wqj indicates the 
eigenfrequency of the phonon in the absence of spin- 
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FIG. 5: (Color online) Dielectric loss of the phonons in 
CdCr2 04. Below Tn one can clearly see the splitting of the 
cubic modes and the appearance of new modes (indicated by 
arrows). 



new modes at 358.7 and 371.3 cm~^. The phonon Ti„(3) 
features a rather clear spHtting into two modes at 463 and 
465.4 cm^^ and a weak new mode appears at 500.5 cm~^. 
On the low-energy side of the ri„(4) mode, at about 
602.3 cm~^, an additional peak becomes visible below 
T/v and the main resonance at 609.4 cm~^ exhibits a 
shoulder at 606.8 cm~^. Therefore, we identify 13 IR 
active phonon modes below T/v for CdCr204. Since from 
the factor group analysis of space group lAi/amd only 
ten IR phonons should be observable, we conclude that 
the true space group should be of lower symmetry than 

In ZnCr2 04 the reported commensurate spin structure 
is accompanied by a tetragonal distortion with a con- 
tracted c-axis and the space group /4m2 (#119) has been 
assignedi^ii^ The corresponding factor-group analysis for 
a 1 X 1 X 1 chemical unit cell for ZnCr204 yields the 
following symmetry distribution of F-point vibrational 
modes: 



phonon coupling at K, Cj is a mode dependent scal- 
ing factor of the anharmonic contributions, and tOav = 
423 cm~^ is the arithmetic average of the IR- and Raman- 
active phonon frequencies at room temperature. 

In these AFM Cr spinels magnetic ordering is accom- 
panied by a structural symmetry reduction. Accordingly, 
a splitting of the degenerate cubic Tiu{j) modes and 
the appearance of new modes is expected. The most 
common structural distortion for spinels occurring be- 
low the Neel temperature is tetragonal with space group 
lAi/amd (#141)ii^i^'i^ with the following irreducible 
representations ^ 

r = 4A2u + 6Eu (IR active) 

-I- 2Aig + 3Big + B2g + AEg (Ramau active) 
+ 2Aiu + A2g + 2Biu + 4S2„ (silent) 

Such a symmetry reduction reportedly takes place in 
CdCr2 04, where the incommensurate spin structure is 
stabilized by a tetragonal distortion with an elongated 
c-axis . ^^1^^ Hence, one expects a splitting of each of the 
four cubic modes into a singlet and a doublet and the 
appearance of two more doublet modes. When zooming 
in on the dielectric loss spectra in the region of the four 
modes in CdCr204 in Fig.El splittings of all modes can be 
detected upon entering into the ordered state similar to 
the previously reported case of polycrystalline samplesi^ 
Note that the effects are very subtle for the Ti„(l) and 
Ti„(4) modes and, therefore, we chose a semi-logarithmic 
plot for these phonons. The Tiu{l) mode develops into 
a broad flat feature indicative of a two-peak structure 
(at 148.6 and 150.4 cm~^) and at about 155 cm~^ a new 
mode emerges. The splitting of the second phonon Tin (2) 
into two modes at 256.2 and 366.1 cm~^, which has been 
described earlier in Ref. 0, is strong compared to the 
other phonons and amounts to approximately 10 cm~^. 
Note that both split modes exhibit a shoulder on the 
high-frequency side which can be interpreted as two more 



r = 8B2 + 11E (IR and Raman active) 

+ 6A1 + 3Bi (Raman active) 

+ 3^2 (silent) 

Note, that in this symmetry the 8-B2 and HE modes 
are both Raman and IR active, because of the lack of in- 
version symmetry. Recently, ZnCr204 was studied using 
synchrotron x-rays and neutron diffractionii^ To deter- 
mine the tetragonal structure in detail, the x-ray inte- 
grated intensity of about 140 different superlattice re- 
flections was measured and fltted within /4m2 symme- 
try. Cr ions were found to occupy six crystallographically 
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FIG. 6: (Color online) Raman spectra of CdCr204 for various 
temperatures. The spectrum at 5 K is shifted for clarity. In 
the magnetically ordered state a sharp resonance appears at 
551 cm~^. Inset: temperature evolution of Raman spectra in 
the region of the T2g{2) and T2g{3) modes {T= 5, 10, 50, 100, 
150, 175, 200, 225, 250, 275, and 290 K from top to bottom). 
Additional features are indicated by arrows and their Raman 
shift is given in cm~^. Spectra below 290 K are shifted for 
clarity. 
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distinct sites: four 8i sites and two 16j sites. Note that 
expanding the tetragonal unit cell by x -\/2 x 2 com- 
pared to the cubic (or simply c-axis contracted tetrago- 
nal) unit cell should lead to the appearance of even more 
extra lines in the Raman and IR spectra. 

A x-ray powder diffraction study was also performed 
for ZnCr204 by Kagomiya et al— and a Rietveld re- 
finement below Tjv was performed for the orthorhom- 
bic Fddd symmetry (#70). For this structure the factor 
group analysis gives: 



'1B2U + 



f 4^23 
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Here one would expect 21 IR-active and 15 Raman- 
active phonons. In principle, the number of observed 
phonons in the IR and Raman data could distinguish be- 
tween the proposed space groups. In the IR experiment, 
however, only the splitting of Tin (2) and a new mode at 
553 cm~^ has been observed;^ and no final conclusion 
can be drawn from IR data alone. Therefore, we will now 
turn to the Raman spectra of CdCr204 and ZnCr204. 



2. Raman scattering 

The Raman spectra for CdCr204 at 290 K and at 5 K, 
below T/v = 7.8 K, are shown in Fig.[6l The five expected 
Raman modes for the cubic structure can be identified 
in the high-temperature data at 290 K. The spectrum 
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FIG. 7: (Color online) Temperature evolution of the four cu- 
bic Raman modes. The dashed lines are fits using Eq. ([5]). 
Solid lines are drawn to guide the eye. 
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(spectra are shifted for clarity). 



consists of three phonon triplets of T2g symmetry (134, 
499, and 600 cm~^), one weak doublet of Eg symmetry 
(343 cm~^), and one singlet of symmetry (647 cm~^). 

The temperature dependence of the eigenfrequencies 
of the four intensive Raman- active modes of CdCr204 
are shown in Fig. [7l With decreasing temperature, the 
phonon eigenfrequencies increase, as usually observed in 
anharmonic crystals. The data for T > 150 K were fit- 
ted using Eq. ^ with uJav = 423 cm~^ (dashed lines in 
Fig. [7]). The 123(1) mode shows a purely anharmonic 
behavior while three phonon modes, T2g{2), T2g{3), and 
Aig, reveal negative deviations for T < 150 K. We ob- 
serve no abrupt shifts of the cubic Raman modes at Tjy. 

Based on its energy and on the fact that it has the 
highest intensity among the other very weak lines in the 
Raman spectra, the mode at about 334 cm~^ is identified 
as the cubic Eg phonon mode. This mode is about four 
times as broad as the T2g and Ag phonon lines. Besides, 
a shoulder is identified on its high energy side, which has 
the same frequency vs. temperature behavior as the main 
line. However, it grows faster in intensity upon cooling. 
Higher order multiphonon scattering is observed in the 
Raman spectra, too. We identified lines with energies 
of approximately 670 and 700 cm~^ (not shown here) as 
overtones of the Eg mode and its satellite. 

In the tetragonal lAi/amd phase at temperatures be- 
low T/v one expects a splitting of the triplet modes into a 
doublet and a singlet and the appearance of two more sin- 
glet modes, overall ten Raman active modes. However, 
additional weak and broad features (indicated by arrows 
in the inset of Fig. [6]) appear in the region of the T2g{2) 
and T2g(3) modes already at temperatures far above T/v. 
Only the feature at 551 cm~^ seems to be directly re- 
lated to the low-temperature phase, because it is still 
absent at 10 K, just above the magnetostructural tran- 
sition. When zooming in on the Raman spectra in the 
region of the 12^(2) mode in Fig. [HI a weak split-off res- 
onance at 493 cm"^ on the low-frequency side is clearly 
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FIG. 9: (Color online) Raman spectra of ZnCr204 at selected 
temperatures measured at parallel polarization configuration. 
Spectra below room temperature are shifted for clarity by a 
constant offset. Arrows indicate new phonon lines at temper- 
atures below Tjv. 



seen upon cooling, similar to the feature at 588 cm~ at 
the low-frequency side of T2g(3). 

At present it is not clear where the additional modes 
observable above T/v stem from, but such effects may 
be attributed to a modification of the selection rules 
similar to observations in resonance Raman scattering 
in semiconductors 4i They can be separated into viola- 
tions of long- wavelength selection rules as purely intrinsic 
bulk effects and additional extrinsic or defect contribu- 
tions. Three main processes and their origin have been 
identified ;"" ' "^ ' (i) q-dependent, intraband Frohlich inter- 
actions, (ii) q- independent, impurity-induced effects, and 
(iii) surface electric field effects at a depletion layer. All 
three of the given mechanisms might be relevant in the 
present case, because comparably large intensities are ob- 
served and defects cannot be ruled out completely. Fur- 
ther studies including a test of polarization selection rules 
on samples with different defect contributions under res- 
onant and nonresonant conditions would be needed to 
further clarify the origin of the additional modes. 

Despite the unknown origin of these additional res- 
onances, the Raman data provides valuable informa- 
tion on the low-temperature symmetry properties of 
CdCr2 04, because the low-temperature eigenfrequencies 
of the T2g{2) and T2g{3) modes are identical with two of 
the weak new modes appearing in the IR data below T/v. 
This strongly suggests that below inversion symmetry 
is lost in agreement with the theoretical predictions by 
Chern and coworkers, whose analysis of the spiral mag- 
netic ground state favored a chiral structure with space 
group IAi22M 

We now discuss the Raman spectra of ZnCr204 
depicted in Fig. [9] for selected temperatures below 
and above the magnetostructural transition at T/v = 
12.5 K. Our measurements at room temperature agree 
with previously published data for single crystals of 
ZnCr204.'^'^ As in CdCr204 we observe the three 
triplets of T2g symmetry (183, 511, and 606 cm~^), the 



weak doublet of the Eg symmetry (457 cm~^), and the 
singlet of Aig symmetry (688 cm~^). 

The temperature dependence of the eigenfrequencies 
of the four intensive Raman-active modes of ZnCr204 
are shown in Fig. [TOl The dashed lines are fits using 
Eq. ([5]). Evidently, all modes reveal deviations from this 
purely anharmonic behavior for T<150 K. The discrep- 
ancies are of the order of 10"'^ and positive for all Raman- 
active modes. The anomalous temperature dependence 
smoothly evolves below 150 K, but significantly becomes 
enhanced just below T/v, where a jump of the eigenfre- 
quencies occurs. We observed also an anomalous increas- 
ing of the phonon line intensities below Tn and an abrupt 
breakdown of their width at Tn- Figures fUTa) and (b) 
show the T dependence of the integrated intensity and 
line-width of the T2g{2) triplet mode, respectively. 

Similar anomalies in the intensity were observed in 
related spinels early on and, subsequently, it has been 
shown that this resonant effect is connected with a 
change of the electronic zone structure caused by mag- 
netic ordering.— For example in CdCr2S4, resonant Ra- 
man scattering is related to a red shift in the absorption 
spectrum caused by the magnetic ordering?^ 

The solid line in Fig. [TTIfb) is a fit of the temperature 
dependence of the phonon line width using 



r(T) = Fo 



exp(/ic wo/SfcsT) — 1 



(6) 



with Fo = 3.46 cm~^ and the eigenfrequency wo = 
516 cm~^. The data can be described very well down 
to the Neel temperature by this expression, which incor- 
porates cubic anharmonicity effects.—"^ 
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FIG. 10: (Color online) Mode energy of the phonons in 
ZnCr2 04 as a function of temperature. The dashed lines are 
fits to the high-temperature results using a simple anharmonic 
term, Eq. 
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FIG. 11: Temperature dependence of the integrated intensity 
(a) and the line width (FWHM) (b) of the Tag (510 cm"^) 
triplet mode. The solid line in (a) is draw in to guide the 
eye and reflects in (b) the expected behavior for ZnCr2 04 
according to Eq. 



Experimentally, besides the shifts of the five cubic 
phonon modes, we observe an appearance of eight new 
phonon modes with frequencies 106.7, 164.8, 257.3, 280, 
348.4, 355, 477.6, and 487.0 cm"^ in the Raman spectra 
of ZnCr204 below Tjy (indicated by arrows in Fig. [9]). 
Unfortunately, the identification of a total of 13 Raman 
and six IR active modes does not allow to conclusively 
discard any of the proposed low-temperature symmetries. 
However, we would like to point out that, in contrast 
to CdCr204, we could not identify any phonon mode 
that appears in both Raman and IR spectra below T/v as 
implied by the assignment of the /4m2 symmetry^iiSiiii 
or the additional orthorhombic lattice distortions within 
space group _F222 (#22) evoked to describe the antiferro- 
magnetic resonance properties of ZnCr204J^ Therefore, 
it seems very likely that at low-temperatures inversion 
symmetry is conserved and IR and Raman modes are 
mutually exclusive. This would single out the assign- 
ment of the orthorhombic Fddd symmetry by Kagomiya 
et al^^ as the possibly true space group in the magneti- 
cally ordered state. 



D. Spin-phonon coupling 

In the Cr spinels spin-phonon coupling is considered to 
be the driving force to reheve magnetic frustration by a 
the magnetostructural phase transition. Already above 
the transition, in the correlated paramagnetic state, the 
spin-phonon coupling shows up in the temperature de- 
pendence of the phonon frequencies. In the following 
we will concentrate on the IR active Ti„(2) modes of 
CdCr204 and ZnCr204, because this mode exhibits a 
clear splitting upon entering into the magnetically or- 
dered state [ Fig. IHc) and (d)]. With respect to the 
behavior expected for the frequency shift due to anhar- 
monicity [solid line in Fig.|4][c)], there is a clear deviation 
which becomes stronger on approaching magnetic order- 
ing. Therefore, it has been attributed to the coupling of 
the lattice to spin fluctuations. Considering a nn Heisen- 
berg spin system, Baltensberger and Helman related the 



frequency shift of lattice vibrations in magnetic crystals 
to the magnetic energy of the system, 

Aw«A(si-Sj), (7) 

by introducing the spin-phonon coupling parameter A, 
which depends on the derivatives of the exchange con- 
stants with respect to the coordinates of the magnetic 
ions 4^4^ The coupling parameter A in the paramagnetic 
state can be determined if both the frequency shift of the 
lattice vibrations and the spin-spin correlation function 
are known. To obtain the spin-spin correlation function 
we follow two independent approaches. 

According to the fluctuation-dissipation theorem the 
magnetic molar susceptibility per spin in the disordered 
state can be written as 

rri.n 

where (s™ • s„) represents the spin-spin correlation func- 
tion between the spin Sm and s„. Restricting the sum to 
nn, the spin-spin correlation function is given by 

_ ksTx{T) .(.-t-l) . 
''^"^ - NAg'fil 3 ' 

i.e., we can directly utilize the experimental susceptibility 
data shown in Fig. [21 

An alternative derivation of the nn spin-spin correla- 
tion function can be obtained from speciflc heat datai^ 
Assuming that the dominant part of the magnetic energy 
is given by the nn Heisenberg interaction, one can extract 
the temperature dependence of the spin-spin correlations 
from the magnetic contribution of the speciflc heat via 

{s^ ■ Sj)c„AT) = const + — i- [ C„(7?)d79, (10) 

where &Na is the number of bonds between adjacent 
magnetic ions per mole. For the analysis of the speciflc 




T (K) T (K) 



FIG. 12: (Color online) Temperature dependence of the spin 
correlations {si ■ Sj) derived according to Eqs. ([Oj and l|10p 
and the eigenfrequency ujo of phonon Tin (2) in a scaling that 
reveals a linear behavior of the two quantities at low temper- 
atures. 
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heat, we used the nn exchange coupling Ji obtained from 
the magnetic susceptibility data derived in Sec. IIII Al 

In Fig. [12] we plot {siSj){T) calculated from the mag- 
netic susceptibility (solid line) and from the magnetic 
contribution of the specific heat (dashed line). In ad- 
dition, the eigenfrequencies ujq of the Ti„(2) mode are 
shown as a function of temperature (solid squares). The 
offsets and scalings of the ordinates were chosen to make 
ujQ coincide with the susceptibility data. Comparing the 
curves for the correlation functions of CdCr204, one can 
see that they show a similar evolution and a good agree- 
ment of their absolute values. At low temperatures both 
exhibit a quasi-Hnear behavior. With increasing tem- 
peratures they increase monotonically and reach zero at 
about room temperature. Up to 125 K the phonon fre- 
quency scales very well with the evolution of {siSj). The 
similarity between the correlation curves can be regarded 
as a justification of the procedure to describe and extract 
the lattice contribution to the specific heat in Sec. IIII Bl 
Since, however, (siSj)c„ suffers from the uncertainty of 
modeling the PDOS and depends on the exact determi- 
nation of the nn exchange coupling Ji , we think that the 
use of {siSj)x leads to a more precise determination of A. 

Assuming that Eq. ^ holds, we directly plot the fre- 
quency shift Auj vs. the spin-spin correlation {siSj)^ and 
{siSj)cm for the measured temperatures above T^r in 
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FIG. 13: (Color online) Eigenfrequency t^o as function of the 
spin-spin correlation functions for CdCr2 04 and ZnCr2 04. 
(a), (c) Luo vs. {siSj}^ obtained from the magnetic suscepti- 
bility data, (b), (d) ujo vs. {siSj)c,„ derived from specific heat 
measurement. 



Fig. [131 In Fig- flST a) one clearly sees a linear behav- 
ior for CdCr204 up to about 100 K. The solid line is 
a Hnear fit to the data [Eq. ([7])] and yields a value for 
the spin-phonon coupling constant in the paramagnetic 
state of Xy^ — 3.3 cm~^. At 125 K the eigenfrequencies 
begin to deviate from the regression curve. The inset 
[Fig.fTsFb)] shows ujq against (siSj)c,„, the spin-spin cor- 
relation function calculated from the magnetic specific 
heat. Below 50 K the data can be fitted linearly giv- 
ing a slope of Ac^ — 3.9 cm~^. Above 50 K it starts 
to deviate until an abrupt breakdown above 150 K oc- 
curs. In the lower panels of Fig. [13] we present our results 
for the ZnCr204 spinel. The plot of ujq versus (siSj)^ 
[Fig. [TSlfcl] shows a linear increase up to about 150 K. 
The slope gives a coupHng constant of A^ — 5.0 cm~^. 
When the frequency shift is plotted against {siSj)c^ we 
attain a value of 4.0 cm~^. This drawing reveals a con- 
vincing linear relation up to about 60 K and then starts 
to segregate from the regression curve [see Fig. flSlfd)]. 

A comparable value A — 6.2 cm~^ has been reported 
for ZnCr204 by Sushkov et alfi using {siSj) derived by 
subtracting the specific heat of non-magnetic ZnGa204 
from ZnCr204. Although we could not follow their ap- 
proach as described above, the agreement with our results 
A^ = 5.0 cm~^ and Ac„ = 4.0, is fairly good. 

Keeping in mind that the spin-phonon coupling leads 
to the magnetostructural phase transition the splitting of 
the Tit, (2) modes allows to obtain further estimates for 
A: The magnetostructural transition can be described by 
the spin-Peierls order parameter Ugp = {Si • 5*2 — 5*1 • S4). 
It accounts for the difference in tetrahedral bonds with 
ferromagnetic or antiferromagnetic coupling and vanishes 
accordingly in the undistorted paramagnetic state. The 
splitting between the two modes can then be expressed 
in terms of n^p (Ref. l43) : 



Auj « A(S'i • S'2 - • 5*4) (11) 

In the case of CdCr204, where the dominant distortion 
is a uniform elongation along the c-axis, the singlet mode 
is lower in energy and the order parameter becomes Ugp — 
—9/4 if one assumes a collinear ground state of the spins. 
For the proposed tetragonal distortion of ZnCr2 04 with 
space group /4m2,^^ Ugp equals 9/2 when assuming an 
oversimplified model of a collinear ground state. 

Via this approach A = 4.4 cm~^ and A — 2.6 cm~^ have 
been extracted previously for CdCr204 and ZnCr204, 
respectively!^ We find that the value for CdCr204 com- 
pares nicely to our results derived in the paramagnetic 
phase. This indicates that the proposed order param- 
eter based on a collinear spin state describes the mag- 
netostructural transition correctly. Assuming that our 
value for A for ZnCr204 is valid to describe also the split- 
ting below T/v, one can estimate Ugp « 2, confirming that 
the assumption of a collinear spin configuration does not 
hold for this compound. 
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IV. SUMMARY 

Using a Quantum Tetrahedral Mean-Field model we 
fitted the susceptibilities of CdCr204 and ZnCr204 and 
obtained the nn and nnn exchange coupling constants 
Ji = 14.7 K, J2 = -4.0 K, and Ji = 33.4 K, J2 = 4.4 K, 
respectively. Moreover, we analyzed the specific heat for 
these systems in comparison to nonmagnetic reference 
compounds. From both the susceptibility and the spe- 
cific heat, we extracted the spin-spin correlation function 
in order to relate it to the shift of the IR phonon Tiu{2) 
which shows the strongest splitting at the Neel tempera- 
ture. We argue that the results derived by using directly 
the magnetic susceptibility appear more reliable, result- 
ing in spin-phonon coupling parameters A = 3.3 cm^^ 
and A = 5.0 cm~^ for the Cd and Zn spinel, respectively. 
Additionally, we show that the observed IR and Raman 
spectra of CdCr204 provide evidence for a loss of inver- 
sion symmetry, thus ruHng out the proposed lAi/amd 



symmetry in the magnetically ordered phase. For the 
low-temperature phase of ZnCr2 04 we depicted that the 
proposed symmetries 74m2 and F222, which imply that 
phonon modes are both Raman and IR active, cannot 
be confirmed by our experiments and, hence, favor the 
orthorhombic Fddd symmetry. 
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